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Abstract 



A semi-discrete Lax pair formed from the differential system and recurrence relation for semi- 
classical orthogonal polynomials, leads to a discrete integrable equation for a specific semi-classical 



orthogonal polynomial weight. The main example we use is a semi-classical Laguerre weight to 
derive a third order difference equation with a corresponding Lax pair. 

1 Introduction 



X 

5^ ' The connection between semi-classical orthogonal polynomials and discrete integrable systems is 



well established. The earliest example of a discrete integrable system in semi-classical orthogonal 
polynomials can be attributed first to Shohat in 1939 [16], then second by Freud [10] in 1976. 
However it wasn't until the 1990's, when the focus within integrable systems shifted from continuous 
to discrete, that Fokas, Its, Kitaev, [6] gave this equation a name; discrete Painleve I, (d-Pi). 

Since then, other examples of discrete Painleve equations have been found through exploring 
the recursive structures of different semi-classical orthogonal polynomial families, including semi- 
classical Hermite [15], semi-classical Laguerre [9] and semi-classical Charlier [17] . 

We define an orthogonal polynomial sequence {fn(z)}^Lo with respect to a weight function 
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w{z) on an interval (a, b) as 

/ P n (z)P m (z)w(z)dz = h n S nm , (1.1) 

J a 

with the corresponding recurrence relation 

zP n (z) = P n+1 + S n P n + RnPn-l (1.2) 

for a monic orthogonal polynomial family P n (z) — z 11 +p nn -iz n ~ 1 +p n n-2Z n ~ 2 + ■ ■ ■ From Bochncr 
[2] we know that if {P„(z)} is a sequence of classical orthogonal polynomials, then P n {x) is a 
solution of the second-order differential equation 

where 4>(z) and ^(z) are fixed polynomials of degree < 2 and < 1 respectively, and A„ is a real 
number depending on the degree of the polynomial solution. As a consequence of this the weights 
of classical orthogonal polynomials satisfy a first order differential equation called the Pearson 
differential equation 

^(<f>(z)w(z))=iP(z)w(z) , (1.4) 
az 

when the degrees of 4> and ip satisfy deg (f> < 2 and deg -0 = 1. However when the deg </> > 2 and\or 
deg -0 > 1 then the weight function produces a class of semi-classical orthogonal polynomials. 

Our approach to semi-classical orthogonal polynomials is to make use of the Laguerre method 
|11| (not to be confused with Laguerre orthogonal polynomials), which derives a pair of first order 
differential equations for a general class of orthogonal polynomials, after the reduction of continued 
fractions. The connection with semi-classical orthogonal polynomials occurs because we associate 
the system with a semi-classical weight function uu(x) of the polynomials with the Pearson equation 
(|1.4|) . For convenience we choose to write the Pearson equation in the following form 

W(z)d z w{z) = V{z)w{z) , (1.5) 

where V(z) = tp — cf>' and W(z) = </>. While our aim and approach is different, the Laguerre 
method has been used to find connections with integrablc systems, including continuous Painlcve 
equations, recently Magnus [14] . found a continuous Painlcve equation of the sixth kind from 
the recurrence coefficients of a semi-classical Jacobi polynomial and Forrester and Witte [3 [8] , 
found a Painlcve equation of the fifth kind, also using the Laguerre method, but one that has been 
extended to include bi-orthogonal polynomials. 

Our work will consider a semi-classical Laguerre weight, similar to that used by [5]. The semi- 
classical Laguerre polynomials have not been as widely explored as the semi-classical Hermite 
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polynomials, nor are they as complex as the semi-classical Jacobi polynomials. Thus. Laguerre 
polynomials are an appropriate choice for finding new discrete integrable systems. 

In section 2 we use the so-called Laguerre method to derive the differential system for semi- 
classical (monic) orthogonal polynomials. We show how the compatibility between the differential 
system and the recurrence relation (jf ,2[) leads to a semi-discrete Lax equation [13j . from which 
discrete integrable systems can be derived for specific semi-classical orthogonal polynomial weights. 
In section 3 we choose the semi-classical Laguerre wei ght l (x) = (x - t) a e-( ax+ i x ) which leads 
to a coupled difference system and a corresponding third order nonlinear difference equation. 



2 The Laguerre Method 

We introduce a moment generating function, the Stieltjes function, 

m = [ ^d( X ) (2.i) 

J z-x 

(Stieltjes transform of the orthogonality measure w(x)) then equations for P n can be summarized 
as 

f(z)P n (z) = Pi 1 \(z) + e n (z), (2.2) 

where P^2i(z) is an associated polynomial to P n {z), with degree n — 1. Although e n (z) is not a 
polynomial, we can define it as 

e n (z) = [ ^^w(x)dx. (2.3) 
J z~x 

The polynomials P l ^2 1 {z) 1 as well as the e n (z), satisfy the same recurrence relations (|1.2[) . but with 
P-i(z) = 0. Additionally we have the following relations between P„, P„ and e n 

PrA% - P n -iP?\ = -K-i (2.4a) 

Pn-ltn ~ P-n^n-l = ~h n -l , ( 2 -4b) 

which can be found using the Christoffcl-Darboux identity. Since both P n (z) and e n (z) satisfy the 
recurrence relation (|1.2p we can give an explicit form of P n (z) and e n (z) defined in terms of the 
recurrence relation's coefficients: 

(n-1 \ n-1 \ 

E sj z^ 1 + E E s i s * R o zTl ~ 2 + ■■■ ( 2 - 5a ) 
3=0 J j=l \k=0 / 



3 



Semi-classical orthogonal polynomials may be defined through a differential difference equation 
of the form 

W(z)d z f(z) = V(z)f(z) + U(z) (2.6) 
which comes from considering W(z)(d z f(z)) and the Pearson equation (11.51) . 

W(z)(d z f(z)) = - / ^^dx = - f * ( J-W(z) W ( X )) dx + [ ^dMx) 



(z — x) 2 J dx 



On the first line we assume the first term reduces to zero because of parameter constraints and 
then we have that: 

where U(z) is a polynomial in z. 

2.1 The fundamental linear system for semi-classical orthogonal poly- 
nomials 



We start with the equation (|2.2jl . differentiate it and multiply by W, so that we can then make use 
of the first order linear differential equation (|2.6[) (with the exception, that for this case we will 
consider the x variable to be dominant). 

wfd x p n + (Vf + u)p n = widvP^ + d^) 

WdxPniP^+e^ + VPniP^+e^ + UPl = W{d x P^\+d x e n )P n (2.7) 

We then go about separating the polynomial expression P„2i and e n so we get the following two 
equivalent expressions, which we denote 0„ 

9„ = W{d x P^\P n - d x P n p£l r ) - VPl - VP n Pi% (2.8a) 
= W{d x P n e n -d x e n P n ) + VP n e n , (2.8b) 

where n is a polynomial bounded by a constant. We try the same method again except this time 
we use /P„_i, which is again differentiated and multiplied by W. 

<W/'. : • /'•'./'• . = d x p(% + fl^c-i (2.9) 
VPn-iiP™! + e n ) + UP n P n .x + Wd x P n ^{P^\ + en) = W(d x Pi% + 9 x e n _i)P n 

(2.10) 



Again we separate the polynomial expression P„2i and e„ to get a second object, which will be 
called fi„: 

0„ = WiPndxPWi-Pjti&Pn-J-VP^pWx-UPnPn-! (2.11a) 



M^(e„<9 x P„-i - P n d x e n -i) + Ve n P n -i 



(2.11b) 



We can express both VL n and 0„ in terms of the recurrence coefficients by substituting the expres- 
sions for P n (|2~5a| and e n (|2.5b[) into Q„ (|2.11b[) and 6„ (|2.8b[) . 



e„ = w{x)k 



+ 



- ^S, (n-l)x- 2 + 



+ 



n + 1 



n + 2 \ Z ^ J / + 3 

3=0 



+V{x) x h„ 



3=0 



\ 1 






X 


j x»+ 2 + 



3=0 

?i - 1 \ 

E*)«"" 1 + 

3=0 / 



(2.12a) 



fin = VK(CC) <^ ft„ 



3=0 



(n-l)x"- 2 -(n-2)[X;S / ) i r''- s 

3=0 



3=0 
3 



fc=0 

n-2 /j-i 



+ ■ 



3 )E E^^-^p n " 4 + 

3=1 \fe=0 



/n-1 \ n-1 

- E $ * n-1 + E E s ^ - ^ ) *"~ 2 + 

\3=0 / 3=1 \fc=0 



FT + ( E S - 

1,3=0 



3=0 



\li 



(2.12b) 

These definitions will be particularly useful when we arc looking at examples of specific semi- 
classical weights. 



^ + (5> )^W+E ^+ i+ E 5 ^ 

\j=0 / 3=0 V fe=0 / 

/n— 2 \ n-2 /j — 1 

\3=0 / 3=1 \fe=0 



Since the recurrence relation (| 1 . 2[) can be expressed in a matrix form 

X S n Rn 



1 







ip n {x), where ip n {x) 



Pn(x) 
Pn-l{x) 



(2.13) 



we collect the important relations we have derived so far and put them in a matrix form so that our 
intended differential system can be written as one expression. We begin with the two expressions 
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(|2.8a|) and (|2.11ap , written in matrix form: 

Pn -P^ ) \ Wd x P n ) \ 6„ + VPnP^ + UPI 

which can easily be solved making use of (|2.4a[) to give: 

Wd x P^\ \_ l I P { n\ -P^2 \( V n + VPn-iP™! + f/P„P„-l 
Wd x P n j h n-l y P„ -P„_j J y 6„ + VP n P£2 X + XJPl 

so that we have two differential equations: 



(2.14a) 



(2.14b) 



W3 x P n = -2—{Q n P n -Q n P n ^), (2.15a) 

wd x p^\ = (si n p£2 l -Q n p£2 2 + vh n -ipM 1 + uh n - 1 p n ). (2.i5b) 

Looking for a second differential relation for P n , we take (|2.15a[) with a reduced index in conjunction 
with the recurrence relation (|1.2p . which leads to 

W{d x P n - X ) = -2— (n„_lP„-l - | ! ^ i ((a; - S„-l)P„-l - P n )) . (2.16) 
hn-2 \ n.n-i J 

However we have no expression to remove the x from the equation, so we consider the problematic 
part of the expression: (x — S n )Q n = (x — S n ) (W(t n d x (P n ) — d x {e n )P n ) + Ve n P n ), which we 
expand using fll.2|) and the differential of (|2.4b[) to get: 

(x - S n )e n = W(-d x e n (P n+1 + RnP„-i) + d x P n (e n+1 + P„e„_i)) 
+VP n (e n+ i + RnCn-x) 
= Q n+1 + R n n n + Vh n (2.17) 

This allows us to remove x from (|2.16|) to give a second differential equation. 

W8 x P n -l = -J—(® n - X P n - n n P n _!) - VP n -l (2.18) 

n>n-i 

We now have a differential system 

1 / U n (x) -0„ (x) 

Wd x i>(x) = Mx) (2.19) 

h n-i \ e n ^(x) -(Q n {x) + V(x)h n ^ 1 ) 

Pn(x) \ 

where tp(x) = . Thus if we give the recurrence and differential equations in a semi- 

Pn-l{x) j 

discrete Lax representation we have 

ip n+ x(x) = L n (x)tp n (x) (2.20a) 
d x ip n (x) = M n {x)i) n {x) (2.20b) 



where 

_ / X — S n ~R n \ 1 | Qn(x) -@n(x) 

\ 1 y W&n-l y n _i(x) -(fln(x) + V(x)ftn_i) 

Here we have identified the Lax matrices I/„ and M„. So given a particular semi-classical weight 
we can identify the polynomials V and W, which in turn lead to expressions for and Q. 



2.2 Compatibility relations 

We now use the differential system (|2 . 19[) with the matrix form of the recurrence relation (|2.13[) in 
order to create a compatibility relation so that relations between Q n and n can be derived. Thus 
we consider the compatibility between the semi-discrete Lax pair, which leads to the semi-discrete 
Lax equation 

d x L n = M n+1 L n - L n M n . (2.21) 

Equating this expression 

1 \ 1 / Q n+ i(x) -Q n+ i(x) \ I x-S n -Rn 




oo/ Wh n y Qn ( x ) -(n n+1 (x) + v(x)h n ) )\i o 

1 I x-S n -Rn \ I Q n (x) -8„(i) 

Whn ~! \ 1 I \ e n _i(a:) -(SlnW + VWh*-! 



we can identify two distinct relations 



(2.22) 



( X - Sn )(^-^-) = R n+1 ^±l-R n ^l + W (2.23a) 
(x-S n )^ = ^± + -^ + V, (2.23b) 

which we can identify as being comparable with the Laguerre- Freud equations pQ . 

Remark 2.1 We should point out that this system could be explored independent of orthogonal 
polynomials by simply setting V = vq+vix+i^x 2 + . . .+v n x n and W = WQ+W1X+W2X 2 + . . .+w n x n 
(where the Vj,Wj are constants) and then see what difference equations are produced for different 
orders of V and W. However since we are interested with the connections with semi-classical 
orthogonal polynomials, we will present a semi-classical weight and then determine V and W. 
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3 A coupled difference equation and corresponding third 
order nonlinear equation 



This method can be demonstrated by using a semi-classical weight Iq(x) = (x — t) a e~( ax+ i x ' 
synonymous with the (associated) Laguerre orthogonal polynomials l(x) = x a e~ x . Our choice of 
deformations for this weight, involve altering the order of the polynomial in the exponential. 

We first consider a deformation in the exponential part of the weight function, the semi-classical 
weight w(x) = (x — t) a e~( ax+! i x ) with a, a, b > and where the support S is an arc from (t — + oo). 
Then from the Pearson equation, we have 

V(x)=a-(a + bx)(x-t) , W(x)=x-t. (3.1) 

and from the consistency relations we have two non-trivial equations 

b(R n+1 + R n ) = -3 n [bS n + (a - bt)] + (2n + 1 + at + a), (3.2a) 
R n+1 [b(S n+1 + S„) + (a - bt)] - R n [b{S n + S n -x) + (a - bt)] = S n - t. 

(3.2b) 

We consider this to be a nonlinear system in terms of the recurrence coefficients R n and S n , which 
has the linear system (|2.20p with the Lax pair: 

(3.3a) 




- bRn (bx + a + b(S n - t))R n 



x - t \ _q )X + a + b(S n -i - t)) bx 2 + x(a - bt) + bR n -n-a-at 



(3.3b) 



for the associated semi-discrete Lax equation (|2.21[) . This system can be called a discrete integrablc 
system due to the existence of the corresponding linear problem, i.e., the Lax pair. In the strictest 
sense we cannot call the equation Painleve since it is third order however there are other examples 
of third order nonlinear difference equations that arc known to be intcgrable [1] . 

Writing this system in matrix form 

b b \ I R n+1 

b(S n+1 + S n ) + (a - bt) -[b(S n + S n -x) + (a-bt)] j V Rn 

-S n (bS n + (a - bt)) + (2n + l + at + a) 

8 



(3.4) 



allows us to solve the system in terms of R n +i and R n and hence we can find a third order difference 
equation in S n 

US n + S„_! + - i)){S n (S n + (f - *))) - J( 2 " + 1 + at + a)((5„ + S n _0 
+ (f - *)) - ( 5 « - *)} x {- 2 (-Sn+i + | - *) - (5« + 2S„+i + 5„ +2 )} 
= |(S*„ +2 + S n+1 + (I - «) )(5„+i(5„+i + (f ~ *))) - J( 2n + 1 + ^ + a)((S'n+2 + 5„+i) 
+ (| - *)) - (S n+1 - t)} x {-2 (S n + | - t) - (5 n _! + 25„ + 5 n+1 )} . (3.5) 

Alternatively by letting 5* n = Q n — Q n -\ in (|3.2p . we are led to an alternative third order 

difference equation in Q n that includes an extra parameter c 

b f Qn - {n + l)t + c Q n -i-nt + c \ 

\a-bt + b(Q n+1 - Q n -i) a—bt + b(Q n - Q„_ 2 ) / 
= -(Qn-Qn-i)[b(Q n -Qn-i)+a-bt] + (2n + l + at + a) , (3.6) 

(where c is the constant of integration). 

4 Conclusion and Outlook 

Given a class of semi-classical orthogonal polynomials (Hermite, Laguerre and Jacobi), we can 
identify a semi-discrete Lax pair and thus an associated discrete integrable system. Using the 
Laguerre weight lo(x) = (x — t) a e'^ ax+ ^ x2 \ we found a new coupled discrete integrable system, 
which is first order in R and second order in S. Combining the two equations gives a third order 
difference equation in S or a third order difference equation in the new variable Q. Since we were 
only interested in the connections with discrete Painleve equations, we have omitted to look at the 
additional t-differential equation (which appears as a consequence of the t parameter in the weight 
function). It is likely that we could use the t-differential equation to find a continuous Painleve 
equation related to our semi-classical Laguerre weight. 

While we have chosen to look at a simple deformation of the classical orthogonal polynomial 
weight associated Laguerre, this scheme can of course be used for deriving a multitude of discrete 
integrable systems, through choosing an appropriate classical weight function. While we have made 
some progress in this regard, finding the corresponding continuous equation is not always possible. 
Working with further examples has shown that when looking for continuum limits, the choice of a 
semi-classical weight function must be of a particular form. 

In this paper we have applied the Laguerre method to a family of classical orthogonal poly- 
nomials, but we expect it is possible that the Laguerre method can be used with other classes 
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of orthogonal polynomials, such as the discrete, multiple or q-orthogonal polynomials. We would 
need to alter the method appropriately, such as choosing an analogue for the Pearson equation 
(since we arc using the "classical Pearson equation" for classical orthogonal polynomials). Thus, 
we could derive a similar scheme for q-orthogonal polynomials given the q-Pearson equation, where 
a natural extension of this would be to consider the q-Laguerre orthogonal polynomials. 
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